The anomalous behavior of coefficient of normal restitution in the oblique impact 
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The coefficient of normal restitution in an oblique impact is theoretically studied. Using a two- 
dimensional lattice model for an elastic disk and an elastic wall, we demonstrate that the coefficient 
of normal restitution can exceed unity and has a peak against the incident angle in our simulation. 
We also explain this behavior based upon a phenomenological theory. 

The coefficient of normal restitution e is introduced to 
determine the normal component of the post-collisional 
velocity in the collision of two materials. The coefficient 
e is denned by 



v(t c 



-ev(0) 
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where v(r) is the relative velocity of the centers of mass 
of two colliding materials at time r measured from the 
initial contact, r c is the duration of a collision, and n 
is the unit vector normal to the contact plane. Though 
some text books of elementary physics state that e is 
a material constant, many experiments and simulations 
show that e decreases with increasing impact velocity Q. 
The dependence of e on the low impact velocity is theo- 
retically treated by the quasi-static theory 

EHEQ- We 

also recognize that e can be less than unity for the normal 
impact without the introduction of any explicit dissipa- 
tion, because the macroscopic inelasticities originate in 
the transfer of the energy from the translational mode to 
the internal modes such as the vibrations 3, 5, 6]. 

While e has been believed to be less than unity in 
most situations, it is recently reported that e can exceed 
unity in oblique impacts 0,]^, In particular, Louge 
and Adams |9| observed oblique impacts of a hard alu- 
minum oxide sphere on a thick elastoplastic polycarbon- 
ate plate and found that e grows monotonically with the 
magnitude of the tangent of the incident angle 7. In 
their experiment, Young's modulus of the plate is 100 
times smaller than that of the aluminum oxide sphere. 
They also suggested that e can exceed unity for the most 
oblique impacts. 

In this letter, we demonstrate that our two- 
dimensional simulation of the oblique impact based on 
Hamilton's equation has yielded an increasing e with 
tan 7 and e exceeds unity at the critical incident angle. 
Finally, we explain our results by our phenomenological 
theory. 

Let us introduce our numerical model [Io|. Our model 
consists of an elastic disk and an elastic wall (Fig. 
The width and the height of the wall are 8R and 2R, re- 
spectively, where R is the mean radius of the undeformed 
disk. The both side ends and the bottom of the wall are 
fixed. We place 800 mass points at random in a disk 
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FIG. 1: The elastic disk and wall consisted of random lattice. 



with the radius R and 4000 mass points at random in a 
wall for the disk and the wall, respectively. We connect 
each mass point with its neighbor mass points by the 
Delaunay triangulation algorithm |ll|. and undeformed 
nonlinear springs are placed on all the connections. 
Each mass point i on the lower half boundary of the 



disk feels the force, F(Z 



'"' — aVo exp(— ali^)!!^ , where 
l ( s '' is the distance between i-th surface mass point of 
the disk and the nearest surface spring of the wall, a = 
300/i?, Vq — amc 2 R/2, m is the mass of each mass point, 

(i) 

c is the one-dimensional velocity of sound, and n s is 
the unit vector normal to the connection between two 
surface mass points of the wall • We should note that 
the strong repulsion F(li 1 ^) is introduced to inhibit the 
penetration of the disk to the surface of the wallQ- Thus, 
the dynamical equation of motion for each mass point i 
of the lower half boundary of the disk is described by 



m 



d 2 r t 



dr 2 



{-ka^ij - fc 6 xf J .}+e(^-ZW)aV exp(-aZW)nW, 

(2) 

where is the position of i-th mass point, Ni is the num- 
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ber of mass points connected to i-th. mass point, x,y is 
the relative deformation of the spring between i-th and 
j-th connected mass points, k a and k b — k a x 10~ 3 /i? 2 
are the spring constants. Here we introduce the step 
function Q(x), i.e. 9(x) = 1 for x > and <d(x) = 
for a; < 0, and the threshold length l t h which is the av- 
erage of the natural lengths of the springs of the disk. 
For internal mass points, the last term of the right hand 
side of eq.(0 is omitted. In most of our simulations, 
we adopt k a = k^ = 1.0 x mc 2 / R 2 for the disk and 
k a = ka = l.Ox \0~ 2 mc 2 / R 2 for the wall. We do not in- 
troduce any dissipative mechanism in this model. Thus, 
during a collision, a part of initial translational energy of 
the disk is distributed into the vibrational energy of the 
disk and the wall. It should be noted that the macro- 
scopic dissipation can be interpreted as the irreversible 
transfer of the energy from the translational motion to 
the internal vibration. When we introduce explicit dis- 
sipations in the model and add the gravity to the disk, 
we have confirmed that the compression of the disk can 
be described by two-dimensional Hertzian contact theory 
0. 

In this model, the roughness of the surfaces is impor- 
tant to make the disk rotate after a collision 1 Oj . We 
modify the positions of the surface mass points of the 
flat wall and the smooth disk by using normal random 
numbers whose average and standard deviation are and 
3 x 10~ 2 i?, respectively. 

Poisson's ratio v and Young's modulus E of this model 
can be evaluated by adding the external force to stretch 
the rectangle of random lattice numerically. We obtain 
Poisson's ratio v = (7.50±0.11) x 10~ 2 and Young's mod- 
ulus E = (9.54 ± 0.231) x 10 3 mc 2 /i? 2 , respectively^. 

We solve the dynamical equation of motion J2J) for each 
mass point with the initial speed |v(0)| = 0.1c and the 
incident angle 7, and determine e for each 7 according 
to eq.Q. All the results in this letter are obtained by 
averaging the results of 100 disks with different configura- 
tions of mass points. We use the fourth order symplectic 
integrator with the time step At = 10~ 3 i?/c. 

Figure is the relation between e and ^1 = — (v(0) • 
t)/(v(0) • n) = tan 7, where t is the unit vector vertical 
to n. The cross point is the mean value and the error bar 
is the standard deviation of 100 samples for each 7. This 
result shows that e increases with increasing ^1 to exceed 
unity, and has a peak around "fi = 6.0. The behavior of 
e having the peak is contrast to that in the experiment 
by Louge and Adams 0- 

Here, let us explain our results. Louge and Adams |9j 
suggest that their results can be explained by the rota- 
tion of the normal unit vector n arising from the local 
deformation of the wall's surface. Thus, we aim to deter- 
mine the angle of rotation of the unit vector a at each 7 
from the theory of elasticity. 

Figure |2| is the schematic figure of a hard disk moving 
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FIG. 2: Numerical and theoretical results of the relation be- 
tween ^1 and e. 
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FIG. 3: The schematic figure of a hard disk sliding on a soft 
wall, x coordinates of both ends of the contact area AB are 
x = x a and x = xt. 



from the left to the right on a wall, where I = \xb — x a \. 
Assuming that I is small compared to R, the perimeter 
of the contact is approximated by 



f{x) = {x-x c ) 2 /2R-y c , 



(3) 



where (x c , y c ) is the lowest position of the disk. To cal- 
culate tana = {f{xb) — f{x a ))/l, we need to know the 
ratio of \x c — x a \ to I. From the theory of elasticity |l3|. 
this ratio can be estimated as 



with 



1 

— arctan — ; -, (4) 



where ^ is the coefficient of the friction. We evaluate \i 
by fi = |J-t|/|J-n| with J = M(v(r) - v(0)), where M is 
the mass of the disk. The cross points in FigQ] represent 
/i calculated from our simulation, where /i has a peak 
around = 3.0. From eqs.|0 and tana is given by 



tana 



26»-l \x r 



2-29 



R 



(5) 
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where we evaluate \x c — x a \ — 0.55-R from the maximum 
value of the compression of the disk. From eqs.Q and 
@, we obtain the relation between and tana. 

Next, let us calculate e from the relation between 
and tana. We introduce the rotated unit vectors, n a 
and t Q , as n Q = cos an — sin at and t Q = sin an + cos at, 
respectively. By introducing e a = — (v(r c ) • n Q )/(v(0) • 
n Q ), we can express e in terms of e a as 



e a + *2 tana 
1 - 4^ tan a ' 



(6) 



where = -(v(0) • t Q )/(v(0) • n a ) and #g = -(v(r c ) 
t Q )/(v(0) • n a ). also can be rewritten as 



(*i - tana)/(l + \&i tana) 



(7) 



On the other hand, in the oblique impact of slipping 
disks, ^2 1S given by 

= - 3(1 + e a )p a (8) 

in the two-dimensional situation^. In eq.ljHJl, de- 
fined by p a = J • t Q |/|J • n a |, is given by 



Ma 



/j + tana 
1 — fi tan a 

To draw the solid line in FigO at first, we calculate p 
and tana for each ^i. By choosing a fitting parame- 
ter e a — 0.95 and substituting eqs. 0-© into eq.JHJ we 
obtain e as a function of a, p, and ^S±. All points are 
interpolated with the cubic spline interpolation method 
to draw the solid curve. Such the theoretical description 
of e is qualitatively consistent with our numerical result, 
though the theoretical value is a little smaller than the 
observed value. 



0.065 ■ 

0.06 
0.055 

0.05 
H 0.045 

0.04 
0.035 

0.03 
0.025,: 




1 2 3 4 5 6 7 



FIG. 4: Numerical and theoretical results of the relation be- 
tween \fq and fi. 

Now, let us consider how p depends on based on 
a phenomcnological argument. For simplicity, we re- 
place the roughness on the surface of the wall by a pe- 
riodic array of asperities. When the disk hits one as- 
perity, a fraction of the energy is absorbed in the wall. 



We assume that the tangential velocity ^ changes 
to vf 1 = (1 — r))v^~ lS> when the disk hits the z-th as- 
perity. Thus, if the disk contacts N asperities during 
the impact, the tangential speed at the release point be- 
comes v t (T c ) — v[ N ^ = (1 — rj) N v t (0). Here we estimate 
the number of contacted asperities during the collision 
as N — pi si, where p is the number of the asperities 
in a unit length on the surface and l s i is the length 
of sliding which can be evaluated as l s i = v t {Q)r c with 
t c = tt(R/c) y/la (4c/v n (0)) Thus, the tangential im- 
pulse Jt = J ■ t is approximated by 



Jt = M{(i- v ) 



phi 



l}v t (0) ~ - V Mpl sl v t (0), (10) 



where we assume small r\. Now, we should answer the 
question how to determine r\. When the impact velocity 
is large enough, r\ becomes smaller because the asperities 
are broken down when the disk hits them. Therefore 
we may assume the form rj = T)o/(l + /3vf (0)/w„(0)) = 
77o/(l + ft'&x) with the introduction of the dimensionless 
parameters 770 and [3. Taking into account J n = J • n — 
— M(l+e)v n (0) and the definition of p = Jt/Jn we obtain 



TrrjopR 

l + e 1 + ^Vl 
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Here we use cos 7 



VW + *i), 7 



(11) 



y/W[/(T+Wf). We adopt the numerical results for the 
value of e at each *i (FigEJl and |v(0)|/c = 0.1. The solid 
curve in Figgis ea. pi|l with fitting parameters j3 = 0.21 
and r/opR = 0.18, which reproduces our numerical result. 
We do not claim that our simple argument explains the 
experimental result because of two fitting parameters. 
However, we would like to emphasize that our picture 
captures the essence of the behavior of p in the oblique 
impact. 

Let us discuss our result. First, we emphasize that 
the novel phenomena of e exceeding unity are obtained 
from the local deformation of the soft wall for the oblique 
impacts of a hard disk. When we simulate the impact 
between a disk and a hard wall, for ki w ^ = 10 x ki d \ 
e fluctuates around a constant to exceed unity abruptly 
around ^1 = 4.5 [l^. This tendency resembles the exper- 
imental results by Calsamiglia et al.^. Thus, for smooth 
increase of e to exceed unity, the wall should be softer 
than the disk. In addition, it is important to fix the ini- 
tial kinetic energy of the disk. So far, we have confirmed 
that e does not exceed u nity when \&i is controlled by 
changing v t with fixed v n |lfl|. 

The second, the initial velocity of the disk and the lo- 
cal deformation of the wall are much larger than those in 
the experimental ones in ref. Q . They cause the most sig- 
nificant difference between our result and their result [9j. 
Because of the high speed impact in our case, there is 
a peak of e for small 7. In fact, our simulation with 
|v(0)| = 0.01c shows the shift of the peak for larger 
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7. Therefore, we expect that our model reproduces the 
result of ref. 9] for the low impact speed. In addition, 
we have carried out simulations with a disk of 400 mass 
points and a wall of 2000 mass points to check the effect 
of the system size. Although there is a slight difference 
between the results, the data are also well reproducible 
by our phenomenological theory. 

The third, the local deformation of the wall also affects 
the relation between n and ^i. In early studies, it has 
been shown that /i depends on the impact velocity [fllllf)l|. 
In our simulation, /x has a peak around \&i = 3.0. This 
behavior is interpreted as that the asperities are flattened 
for large Vt- Equation (|f 1 1 indicates that /i can increase 
with increasing 7 as in ref. j9j if we choose a suitable set 
of j3 and 7/0 p. The difference between the results of their 
experiment and our simulation may be explained by the 
choice of these parameters. 

The fourth, we adopt the static theory of elasticity to 
explain our numerical results in this letter. However, it 
is important to solve the time-dependent equation of the 
deformation of the wall's surface to analyze the dynamics 
of impact phenomena. The dynamical analysis will be 
our future task. 

In final, we indicate that the friction coefficient is de- 
rived from our Hamiltonian model. This friction comes 
from the irreversible energy transfer of the macroscopic 
translational motion to the internal vibration. The irre- 
versibility is indeed related to the second law of thermo- 
dynamics if the number of the internal degrees of free- 
dom is infinite. However, the irreversibility in our sys- 
tem which includes only 10 3 mass points in a disk is in- 
complete as we can see that the Hertzian contact theory 
cannot be recovered without introduction of the explicit 
dissipation. We believe that the complete treatment of 
inelastic collision of macroscopic materials without intro- 
duction of dissipations will be a fundamental subject of 
nonequilibrium statistical mechanics. 

In summary, we have carried out the two-dimensional 
simulation of the oblique impact of an elastic disk on an 
elastic wall. We have found that e can exceed unity in 
the oblique impact, which is attributed to the local de- 
formation of the wall. We have estimated the magnitude 
of the local defomation a based on the static theory of 
elasticity and derived the relation between e and ^1 by 
taking into account the rotation of the normal unit vec- 
tor of the wall's surface. The relation between and \&i 



is also related to the local deformation and is explained 
by the simple argument. 
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